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The gluon and ghost propagator and the influence of Gribov copies * 
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The dependence of the Landau gauge gluon and ghost propagators on the choice of Gribov copies is studied in 
pure SU{3) lattice gauge theory. Whereas the influence on the gluon propagator is small, the ghost propagator 
becomes clearly affected by the copies in the infrared region. We compare our data with the infrared exponents 
predicted by the Dyson-Schwinger equation approach. 



The non-perturbative behaviour of the gluon 
and ghost propagators in Yang-Mills theories is of 
interest for the understanding of the mechanism 
of confinement of gluons and quarks. In particu- 
lar the infrared behaviour of the ghost propaga- 
tor in the Landau gauge is related to the so-called 
Kugo-Ojima confinement criterion 0, which ex- 
presses the absence of coloured massless asymp- 
totic states from the spectrum of physical states 
in terms of the ghost propagator at vanishing mo- 
mentum. On the other hand the suppression of 
the gluon propagator in the infrared was argued 
to be related to the gluon confinement 0|- 

Zwanziger |3| has suggested that the behaviour 
of both propagators in Landau gauge results from 
the restriction of the gauge fields to the Gribov 
region, where the Faddeev-Popov operator is non- 
negative. Generically, one gauge orbit has more 
than one intersection (Gribov copies) within the 
Gribov region. In this contribution we assess the 
importance of this ambiguity for the ghost and 
gluon propagator in SU (3) gauge theory on a fi- 
nite lattice. 

In the continuum the gluon and ghost propa- 
gators have been computed from a coupled and 
truncated set of Dyson-Schwinger (DS) equa- 
tions 0] . Representing the gluon and ghost prop- 
agators in the Landau gauge as 
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the low-momentum behaviour of the correspond- 
ing dressing functions was found to be closely 
constrained by the respective infrared exponents, 
Zgi oc (g^)^'", Zgh oc (q^y with 0.5 < K < 1. 

There have been only a few numerical lattice in- 
vestigations of the SU{i) ghost propagator in the 
past |3| contrary to the gluon propagator. Here 
we present data for both propagators in Landau 
gauge measured at the same ensemble of SU(i) 
gauge- fixed configurations and concentrate on the 
infrared behaviour inside the Gribov region. The 
strategy of the investigation is similar to a previ- 
ous analysis of the SU{2) ghost propagator 0]. 

The Landau gauge condition is implemented 
by searching for a gauge transformation ^Ux.fi = 
gxUx,fj,gl_^^f^ which maximises the gauge func- 
tional Fu [g] oc ^ 9^e Tr ^Ux,fj. where the gauge 
field Ux,fj. is provided by Monte Carlo simula- 
tions. This functional has many local extrema 
whose number increases with the lattice size and 
decreasing f3 = G/g"^. All Gribov copies {^U} be- 
long to the gauge orbit created by U and satisfy 
the lattice Landau gauge condition dfj^^A^^^ = 
with 
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Before investigating the infrared behaviour of 
both propagators we first check the influence of 
different ways to select Gribov copies on the av- 
erage propagators taken at some low momenta. 

As it was previously found the Gribov copy 
ambiguity affects the measurements of the ghost 
propagator in the case of SU (2) 0,0 ■ Also from a 
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Figure 1. The ghost and the gluon propagator 
Z{q^)/q^ in lattice units at the two lowest momenta 
q{k) as a function of (3 for the first (fc) and best (be) 
gauge copies (lattice size 24*). 

recent study the SU{?>) gluon propagator is seen 
to depend systematically on the choice of Gri- 
bov copies ;Sj contrary to common belief. We 
report on first measurements on 16* and 24* lat- 
tices at /3 = 5.8, 6.0, and 6.2 as well as 32* at 
(3 — 5.8. In the update the usual heatbath algo- 
rithm with microcanonical overrelaxation steps is 
used. Then for each configuration Ncp — 30 ran- 
dom copies on a 16* lattice {Ncp — 40 on 24* 
and Ncp = 10 on 32*) have been gauge-fixed 
using the standard overrelaxation method until 
maxa;(9^?42;^^)^ < 10^^* has been reached. 

The propagators were evaluated on the first 
(fc) and the best (be) (with respect to the 
gauge functional) gauge- fixed copies. The gluon 
propagator was calculated for all momenta us- 
ing the fast Fourier transformation of the gauge 
potentials The ghost propagator was ob- 

tained using the conjugate gradient method 
with a plane wave source only for the lowest 
non-vanishing momenta. For each momentum 
squared, g^(fc) = (4/a^) Yl,^ sin^ (irk^/L^), the fi- 
nal propagator is an average over the respective 
ensemble of gauge-fixed copies (fc) or (be) and 
over different k giving rise to the same q^{k). Fol- 
lowing Ref. only momenta with k satisfying 
T^p, ~ 1/4 • (Yip ^t^Y ^ 1 have been used in the 
analysis of the gluon propagator. 

The results from a 24* lattice at the chosen val- 
ues of f3 are shown in Fig. ^ for some lowest mo- 



menta. The ghost propagator is clearly affected 
by the Gribov copy problem. The effect is big- 
ger than the statistical error and increases with 
decreasing momenta. On the other hand, for the 
lowest momentum of the gluon propagator the 
impact of Gribov copies is completely below the 
statistical error. For other momenta there are 
some effects visible as shown in the lower half of 
Fig.H 

The upper half of Fig. [21 shows both propaga- 
tors for the best copy as a function of the momen- 
tum q scaled with an appropriate lattice spacin; 
a. In order to compare to other studies 
we have used a~* = 1.53, 1.885 and 2.637 GeV 
for (3 = 5.8, 6.0 and 6.2, respectively. The lower 
half shows the ratio between the mean value on 
the first and best copy of the respective propa- 
gators dressing function, {Z'-^'^'>)/{Z^^'^^). Due to 
the low number Ncp of gauge copies such ratios on 
a 32* lattice are not shown there. For the lowest 
momenta the ghost propagator is systematically 
overestimated if measured on an arbitrary (first) 
gauge copy, while for the gluon propagator the 
statistical noise is dominant. However, so far we 
have only measured the ghost propagator for a 
single k without permutations of components. 

As in ,6| we also found outliers in the data 
of the ghost propagator at the lowest momenta 
k — (1, 0, 0, 0), shown in Fig. |31 which may make 
statistical analyses difficult, if arbitrary (first) 
gauge copies are studied only. However, the rea- 
son is unknown yet. 

To parametrise the infrared behaviour, we fit- 
ted the measured dressing function in the form 
Zgi = Cgi(g2)2« and Zgh = Cgh{q^y with com- 
mon K to the data at some low q values. For the 
best fit with (/max = 0.43 we found n = 0.23(1) 
for the best copy, far away from the DS predic- 
tion. However, despite k tends to increase as gmax 
decrease a reasonable statement cannot be made 
from the lattice sizes L* used. 

In summary, we have reported first results 
studying the influence of the SU{i) ghost and 
gluon propagators on Gribov copies. We have 
found that the ghost propagator is more affected 
than the gluon propagator. Measuring the ghost 
propagator on an arbitrary first gauge-fixed con- 
figuration the propagator is systematically over- 
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Figure 2. The upper figures show the dressing functions of the ghost and gluon propagator measured on the best 
copies as functions of the momentum q scaled to physical units for /3 = 5.8, 6.0 and 6.2 and various lattice sizes. 
The lower ones show the ratio determined from the first (fc) and best (be) gauge copies. 



estimated, the effect is largest for the lowest mo- 
menta. On the other side there is an irregular im- 
pact on the gluon propagator. Concerning the in- 
frared behaviour the results from the lattice sizes 
used so far do not allow to confirm the proposed 
scaling behaviour using the ansatz discussed with 
an infrared exponent k > 0.5 simultaneously de- 
scribing both propagators. 
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Figure 3. The ghost propagator at fc = (1,0,0,0) 
vs. the first (fc) and best (be) gauge copy on a 24* 
lattice. 
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